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Abstract 

We report a combination of the CDF and DO measurements of the inclusive single 
top quark production cross section in the s- and i-channels, <7 s +t, in pp collisions at a 
center of mass energy of 1.96 TeV. The total integrated luminosity included in CDF's 
analysis is 3.2 fb _1 and DO's analysis has 2.3 fb^ 1 . A Bayesian analysis is used to extract 
the cross section from the distributions of multivariate discriminants provided by the 
collaborations. For a top quark mass mt = 170 GeV/c 2 , we measure a cross section of 
2.76ia47 pb. We extract the CKM matrix element \V tb \ = 0.88 ± 0.07 with a 95% C.L. 
lower limit of \Vtb\ > 0.77. 



1 The Tevatron Electroweak Working Group can be contacted at tev-ewwg@fnal.gov. 
More information can be found at [Ettp777tevewwg.fnal.gov 
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1 Introduction 



In the standard model (SM), top quarks are expected to be produced singly in pp collisions 
through s-channel or t-channel exchange of a virtual W boson [lj. The reasons for studying sin- 
gle top quarks are compelling: the production cross section is directly proportional to the square 
of the CKM matrix [2] element \V t b\, and the cross section measurement is sensitive to fourth- 
generation models, models with flavor-changing neutral currents, and other new phenomena [3]. 
The main backgrounds to the single top final state are jy+jets production, ti production, and 
QCD multijet production. Electroweak production of single top quarks is a difficult process to 
measure because the expected production cross section (a*+ e ° ry ~ 2.9 — 3.5 pb [HE]) is much 
smaller than those of the background processes, and it is also smaller than the uncertainty on 
the total background. The presence of only one top quark in the event provides fewer features 
to use in separating the signal from background, compared with measurements of top pair 
production (ti), which was first observed in 1995 [6]. 



The CDF and DO collaborations have developed a variety of sophisticated multivariate 
techniques that take advantage of the differences in the kinematic distributions between signal 
and backgrounds. These have successfully been used to observe the process, at the 5.0 standard 
deviation level of significance, separately by DO [7] and by CDF [8]. Further details of these 
analyses are given in earlier publications presenting evidence for single top production in smaller 
samples of data QUI EE]. 

Because the presence of single top production has been firmly established by the two col- 
laborations and due to the computational difficulty of computing the significance of the two in 
combination, we do not perform a hypothesis test of the combined results. The cross section 
measurements, however, remain limited in precision by the sizes of the data samples collected 
by the two collaborations, and thus a combination of the results provides a more precise mea- 
surement and also stronger constraints on the CKM matrix element \Vtb\- 



The CDF and DO collaborations both measure the single top cross section and extract |V#,| 
using a Bayesian statistical analysis. We use the same method here to combine the two mea- 
surements. The CDF and DO collaborations use independently developed computer programs 
for performing the calculation, and we have verified that both codes reproduce the individual 
experiment results. Since the two codes produce consistent results we use the DO code here to 
obtain the combined Tevatron results. CDF and DO have exchanged histograms of their final 
multivariate discriminants for the observed data and all signals and backgrounds in each analy- 
sis channel. Furthermore, all systematic uncertainties - those on the predicted rates, histogram 
shapes, and bin-by-bin independent Monte Carlo statistical uncertainties, have been exchanged 
as well. 



2 Method 
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Rate and shape uncertainties are categorized by their source for use in evaluating correla- 
tions, between signals and/or backgrounds within a channel, between channels, and between 
experiments. Sources of systematic uncertainty that are common to the measurements of the 
two collaborations have been assigned 100% correlation, and other sources of uncertainty are 
assumed to be uncorrelated. A list of systematic uncertainties is given in Tableland discussed 
in more detail below. 

The method used for extraction of the cross section chosen by the collaborations is a 
Bayesian technique [12J. The posterior probability density is given by 

p(a s+t ) = J L(dataK +t , {6})ir(a s+t )n({6})d{6} , (1) 

where L is the joint likelihood function of all bins in all histograms 

(V _|_ h.\n ip -{s l +b i ) 

l= n {±± ^ — , (2) 

bins, channels 1 

{9} is the set of nuisance parameters (such as integrated luminosity, fo-tag efficiency, and other 
rate and shape parameters), and 7r({0}) is the product of the priors encoding the systematic 
uncertainties on {6}. The predictions for the number of signal events (sj) and background 
events (fej) depend on the values of the nuisance parameters. 

The prior density for the signal cross section is taken to be a nonnegative flat prior, 7i(<J s+t ) = 
l/cmax for a > 0, and n(a) = otherwise. Here, <r max is the highest signal cross section axis 
point (7 pb for the combination). The priors on {6} are truncated Gaussians, where the 
truncations are introduced so as not to allow negative predictions of s« or hi for any template 
in any bin. The integral over d{8} is performed through numerical integration. 

We quote the measured cross section as the value that maximizes the posterior likelihood, 
and the uncertainty as the shortest interval containing 68% of the posterior area. 

In addition to the above, the DO code uses a Gamma prior for the MC statistics uncertainty, 
reflecting the Poisson nature of MC statistics, and having the additional advantage that an 
analytic solution to the integral in Eq. [1] exists. This does not bias the cross section but 
reduces the computing time considerably. 

The combination procedure includes all channels from both collaborations in the Bayesian 
analysis, including correlations of systematic uncertainties between experiments. We combine 
the CDF super-discriminant, the CDF missing transverse energy (Et) plus jets analysis, and 
the DO combination discriminant. While there are some differences in the detailed treatments 
of uncertainties, the CDF and DO codes give results consistent with each other when run over 
the individual inputs. We have verified this by running the DO code over the CDF inputs. 

Tabled] (Table [2]) lists the measured cross sections with and without systematic uncertainties 
for the two ingredients to the CDF combination as well as the CDF combination result, both 
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Table 1: Results summary for CDF's £+jets super-discriminant analysis, ^T+jets analysis, and 

CDF's total combination, comparing the CDF code to the DO code. This is for the published 
CDF measurements at m t = 175 GeV/c 2 . 

Analysis Cross Section (pb) Cross Section (pb) 

CDF Calculation DO Recalculation 

No systematic uncertainties 

Super-discriminant 2.29±^tl 2.291^1 

^T+jets 4.181{;1 4.191$ 

Combined 2.49l£| 2.481^2 

Including systematic uncertainties 

Super-discriminant 2.llgj 2.1±gj 

^T+jets 4.9l|| 5.211;? 

Combined 2.3±g| 2.3±g-| 



Table 2: Results summary for CDF's £+jets super-discriminant analysis, ^T+jets analysis, and 
CDF's total combination, comparing the CDF code to the DO code. This is for the re-calculation 
corresponding to m t = 170 GeV/c 2 and includes all systematic uncertainties. 

Analysis Cross Section (pb) Cross Section (pb) 

CDF Calculation DO Recalculation 

Super-discriminant 2.171^ 2.12l&j£ 

$ T +jets 5.01U 5.4l|j 

Combined 2.35±g;|g 2.29±g;|2 



for the original CDF calculation and the recalculation using the DO code, for a top quark mass 
of m t = 175 GeV/c 2 {m t = 170 GeV/c 2 ). The DO code reproduces the CDF results, indicating 
that all inputs and uncertainties are being read properly. 

Different conventions have been used by the CDF and DO analyses for the top quark mass 
and the theoretical single top cross section. CDF assumes m t = 175 GeV/c 2 in its extraction of 
the central value of the single top production cross section, while DO assumes m t = 170 GeV/c 2 . 
The current world average top mass is 173.1 ± 1.3 GeV/c 2 [13J. CDF includes in its p-value 
calculation and the |Vtb| limit the uncertainty in m t and thus all rate and shape uncertainties 
have been evaluated for the single top signal templates and the ti background templates. DO 
does not have single top or ti samples for m t = 175 GeV/ c 2 , so CDF has shifted its predictions to 
mt = 170 GeV/c 2 , and thus the results quoted here assume mt — 170 GeV/c 2 . The discriminant 
shapes and normalizations have been shifted for s-channel and t-channel signal as well as for 
ti background. The ti background cross section at m t = 170 GeV/c 2 is 7.91 ± 0.95 pb [14J. 
Other background normalizations that weakly depend on ti have not been modified. This has 
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a small impact on the measured cross section as can be seen in Table [2J 

The other difference in convention is the choice of theoretical cross section to compare the 
measurement with. CDF chooses Harris and Sullivan's next-to-leading order (NLO) calculation 
with cross sections of a* heory = 0.99±0.07 pb and a t theory = 2.15±0.24 pb at m t = 170 GeV/c 2 @], 
while DO chooses Kidonakis' NLO plus soft-gluon corrections calculation with cross sections 
of (T* heor y = 1.12 ± 0.05 pb and af ieovy = 2.34 ± 0.13 pb [5J. This choice has no effect on the 
measured cross section, but it does affect the extraction of \V t b\ and its limits. We report both 
extractions of \Va,\ here. 



3 Systematic Uncertainties 

The complete set of systematic uncertainties and their correlations from the CDF and DO 
analyses have been included in the combination. They are categorized in Table [3j 

The systematic error categories in Table [3] mainly follow those used for the Tevatron top 
mass combination [13J aiming to lump together sources of systematic uncertainty that share 
the same or similar origin as well as the same correlation between the two experiments. Each 
error category is briefly described below. 

Luminosity from detector: That part of the luminosity uncertainty that comes from the un- 
certainty on the luminosity detector acceptance and efficiency. It is uncorrelated between 
CDF and DO. 

Luminosity from cross section: That part of the luminosity uncertainty that comes from 
the uncertainty of the inelastic and diffractive cross sections. It is correlated between 
CDF and DO. 

Signal modeling: The systematic uncertainty arising from uncertainties in the modeling of 
the single top signal. This includes uncertainties from variations in the ISR, FSR, and 
PDF descriptions. It also includes difference in the hadronization models. It is correlated 
between CDF and DO. 

Background from MC: The systematic uncertainty arising from uncertainties in modeling 
of the different background sources that are correlated between CDF and DO. It includes 
ti and diboson normalization uncertainties obtained from theoretical calculations. 

Background from data: The systematic uncertainty arising from uncertainties in modeling 
of the different background sources that are obtained using data-driven methods and 
are uncorrelated between CDF and DO. It includes the uncertainty on the IV+jets and 
Wbb/Wcc normalization, scale factor and shape as well as the uncertainty on the multijet 
modeling and normalization. It also includes the uncertainty due to MC statistics. 
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Detector modeling: The systematic uncertainty arising from the uncertainty on the event 
detection efficiencies for object identification and MC mismodeling of data. It is uncor- 
rected between CDF and DO. 

b-tagging: The systematic uncertainty coming from the uncertainty on the b-tagging and 
mistag rate and shape modeling. It is uncorrelated between CDF and DO. 

dJES: The part of the JES uncertainty which originates from limitations in the calibration 
data samples used. For CDF this corresponds to uncertainties associated with the in- 
dependent JES corrections which are estimated using di-jet data events. For DO this 
includes uncertainties in the calorimeter response for light jets, uncertainties from 77- and 
PT-dependent JES corrections, and other small contributions. It is uncorrelated between 
CDF and DO. 



Table 3: Sources of systematic uncertainty in the CDF and DO single top analyses, together 
with the rate range for each uncertainty and whether the shape effect on the final discriminant 
is considered for each uncertainty or not. 



Systematic Uncertainty 


CDF 


DO 




Correlated between 
the two experiments 




Rate Shape 


Rate 


Shape 




Luminosity from detector 


4.5% 


4.6% 






Luminosity from cross section 


4.0% 


4.0% 




• 


Signal modeling 


2.2-19.5% • 


3.5-13.6% 




• 


Background from MC 


12.1-12.4% • 


15.1 % 




• 


Background from data 


17-40% • 


13.7-54% 


• 




Detector modeling 


0-9% 


7.1 % 






b-tagging 


0-29% • 


2-30% 


• 




dJES 


0-16% • 


0.1-13.1% 


• 





We do not include the small uncertainty due to the top quark mass as our result is computed 
at m t = 170 GeV/c 2 . We also do not include the systematic uncertainty on the single top cross 
section theory prediction when measuring the cross section, but we do include it as an additional 
systematic when extracting \V t b\- 



6 



4 Results 



Figure [T] shows the posterior probability distribution of the combined CDF and DO analyses, 
and Table 0] lists the measured cross section. 

Table 4: Measured single top cross section and \V t b\ value and lower limit for two different 
theory cross sections. 



2.761&I? pb 
0.88 ±0.07 

3.46 pb 0.77 

0.91 ±0.08 

3.14 pb 0.79 



measured cross section 
\V tb \, cr*+ t ory = 3.46 pb 
\V tb \ 95% C.L. limit, (T*+7 y 
|^|,ajT ry = 3.14 pb 
\V tb \ 95% C.L. limit, aT+t iy 




single top cross section [pb] 



Figure 1: The posterior probability distribution of the combined CDF and DO analyses. 



The most probable value of the combined s-channel and t-channel cross sections is 2.76_ ; 47 pb 
for a top quark mass of 170 GeV/c 2 . This value lies in between the CDF measurement of 
2.35^50 pb and the DO measurement of 3.94 ±0.88 pb as expected. The combination improves 
the cross section uncertainty from about 22% for each experiment to 19% for the combination. 
We have also computed the compatibility of the two measurements with each other, based on 
the uncorrelated systematic uncertainties, and find that they are 1.6 standard deviations apart. 
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In order to extract information about |V#,| from the combined data, we assume that the off- 
diagonal CKM matrix elements \Vt s \ and \ Vtd\ are much smaller than |V#,|, but we do not make 
assumptions about the unitarity of the 3x3 CKM matrix, allowing for a fourth generation of 
quarks. We obtain the \ Vtb\ 2 posterior by dividing the measured cross section by the theoretical 
single top cross section for V th = 1 and a top quark mass of 170 GeV/c 2 . 




Figure 2: The posterior probability distribution of the combined CDF and DO analyses for \ Vtb\- 
The right plot shows only the posterior, Vtb measurement and limit using a theory cross section 
of 3.46 pb. 



We obtain \V t b\ = 0.91 ± 0.08 using the theoretical calculation of cr* + e ° ry = 3.14 x \V t b\ 2 pb 
of Harris and Sullivan, and \V t b\ = 0.88 ± 0.07 using the theoretical calculation cr*+ e ° ry = 
3.46 x \V t b\ 2 pb of Kidonakis. We limit \V t b\ > 0.79 at the 95% C.L. assuming a flat prior in 
\V t b\ 2 from to 1 using the calculation of Harris and Sullivan, and limit |V#>| > 0.77 at the 95% 
C.L. using the calculation of Kidonakis. The two \V t b\ posteriors are shown in Figure [2j In 
comparison, the CDF extraction of \V t b\ at 175 GeV/c 2 is 0.91 ± 0.14 and the DO extraction of 
| V^b | at 170 GeV/c 2 is 1.07 ±0.12. The relative \V t b\ uncertainty improves from 14% (CDF) and 
11% (DO) to 8% for the combination, more than for the cross section due to the square root 
dependence. 



5 Summary 

We combine the CDF and DO measurements of the single top production cross section a s+t 
including the effects of all correlated and uncorrelated systematic uncertainties. The combined 
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measurement is cr s+t = 2.76_ ; 47 pb, in agreement with the SM expectation. A summary of the 
cross section measurements is shown in Figure [3j From this we obtain measurements of \Vtb\ 
using two different theory cross section normalizations for a top quark mass of 170 GeV/c 2 . We 
find | Vtb\ = 0.88 ± 0.07 for a cross section of a*+ e ° ry = 3.46 x \V t b\ 2 pb and set a 95% C.L. lower 
limit of \V tb \> 0.77. 



Single Top Quark Cross Section 



CDF Lepton+jets 3.2 fb ' 
CDF MET+jets 2.1 fb"' 
D0 Lepton+jets 2.3 fb 1 
Tevatron Combination 

Preliminary 

H B.W. Harris et al., PRD 66, 054024 (2002) 
H N. Kidonakis, PRD 74, 1 1 401 2 (2006) 



August 2009 



2.17 



+0.56 
-0.55 



pb 



H 5.0 !|| pb 
3-94 li| pb 
2.76 Pb 



'top 



170 GeV 



2 4 6 8 

a (pp -» tb+X, tqb+X) [pb] 



Figure 3: Tevatron single top cross section measurements and their combination. 
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